Binder Parameter of a Heisenberg Spin-Glass Model in Four Dimensions 
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We studied the phase transition of the ±J Heisenberg model with and without a random 
anisotropy on four dimensional lattice LxLxLx(L + l) {L < 9). We showed that the Binder 
parameters g(L, T)'s for different sizes do not cross even when the anisotropy is present. On the con- 
trary, when a strong anisotropy exists, g(L,T) exhibits a steep negative dip near the spin-glass phase 
transition temperature Tsg similarly to the p— state infinite-range Potts glass model with p > 3, 
in which the one-step replica-symmetry-breaking (RSB) occurs. We speculated that a one-step 
RSB-like state occurs below Tsg, which breaks the usual crossing behavior of g(L,T). 
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Recently, the low temperature phase of the XY and 
Heisenberg spin-glass (SG) models has been attacked 
a great interest. It is known that vector SG models 
have a chiral symmetry in addition to the continuous 
symmetry^. Consequently, these models may have both 
SG order and chiral glass (CG) order. Kawamura and 
coworkers claimed that decoupling.^ the spin and chiral 
valuables occurs at long distancesoa, and that the CG 
order is realized at low temperataii(jS|in three dimensions 
(d = 3), but the SG order is notffl'Btl Although the ex- 
istence of the CG order has been accepted, controversy 
exists on the SG order. Maucourt and GrempelEl studied 
the domain wall energy W(L) of the ±J XY model on 
finite d — 3 lattices of L 3 at absolute zero temperature 
(T = 0) and speculated that W(L) may increase with 
the linear size L. The same speculation was also given 
by Kosterlitz and Akinotl. The present authors exam=. 
ined the stiflpess of the ± J Heisenberg model at T = CU 
and T ^ 00, and suggested that the stiffness exponent 
9 in d = 3 has a positive value for T/J < 0.19. They 
also performed a Monte-Carlo (MC) simulation of the 
same model and found that the SG susceptibility xsg 
exhibits a dwergent behavior toward the temperature of 
T/J ~ 0.1 8t3. These facts strongly suggests that the 
SG order also occurs in the XY and Heisenberg mod- 
els in d = 3. This view of the SG order was supported 
in recent studies of the aging-affect of the spin autocor- 
relation function (£L(0)Si(t)p3 and the non-equilibrium 
relaxation of xsgc-i However, the Binder parameters 
g(L, T)'s for different L neither cross nor come together 
in both the XY0 and HeisenbergQ SG models. 

This fact raises an objection against the occurrence of 
the SG order, because it is believed that the most reliable 
evidence of the occurrence of the SG phase transition 
is the . -crossing of g(L, T)'s at the same temperature of 
TsGC-Jiij'tJ- In fact, it was reported that the crossing of 
g(L, T)'s occurs in both the XYt§ and Heisenbergca SG 
models in four dimensions (d — 4). Thus we are faced 
by a serious problem that we have gotten the opposite 
views of the SG order from the studies of the different 



quantities. 

The problem may arise from a poor knowledge of the 
property of g(L,T) of the vector SG model. In the non- 
frustrated system, the crossing of g(L,T)'s occurs at Tq 
with some positive value of g (= g(L, 2c)). This property 
results from the fact that, in the thermodynamic limit, 
5(00, T) = for T > T c and g(oo,T) = 1 for T < T c . 
However, the same will not always be true in the SG 
model, because g(oo,T) for T < Tsg will take differ- 
ent values due to the occurrence of the replica symmetry 
breaking (RBSl_.In fact, it was revealed by Hukushima 



T SG ) of the infinite-range p- 



and Kawamurao that 

state Potts glass modelEi takes different values depend- 
ing on the state number p, where T~ = \im e ^o(T — e). 
Therefore it is crucially important to reveal the property 
of g(L, T) in such vector SG models in which the SG 
phase transition occurs at Tsg 7^ 0. 

In this paper, we report that g(L, T) of the Heisenberg 
SG model in d = 4 exhibits a behavior quite different 
from that of the Ising SG model. We reexamined g(L, T) 
of the ±J Heisenberg models in d = 4 with and with- 
out a random anisotropy of magnitude D. In the case 
of D = 0, when L is increased, the increment of g{L,T) 
with decreasing temperature becomes less steep at low 
temperatures, and the crossing of g(i,T)'s which was 
suggested for small IxB is released. When D ^ 0, g(L,T) 
for each L decreases, particularly around Tsg, a- n d this 
property is enhanced when D is increased, where Tsg is 
the SG transition temperature estimated from the scaling 
plot of XSG- In particular, g(L,T) for large D exhibits a 
negative dip near Tsg which deepens with increasing L. 
These facts are quite interestingj-because the anisotropy 
would stabilize the SG orderBEj. We believe, hence, 
that the absence of the crossing of g(L,T) 's for finite L 
says nothing about the presence of the SG order in this 
model. We speculate that, even for D = 0, g(oo,Tg G ) 
takes a small positive value (or a negative value) due to 
the occurrence of a one-step RSB-like state. We hope 
our findings will help to understand the low temperature 
phase of the Heisenberg SG model in d = 3. 
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We start with the anisotropic ± J Heisenberg model on 
a hyper cubic lattice of L x L x L x (L + 1)(= N) with 
skew boundary conditions along three L directions and a 
periodic boundary condition along the (L + 1) direction. 
The Hamiltonian is described by 



H 



{ij 



i j Si S j 



E 



(1) 



where Si is the Heisenberg spin of | | = 1 and Sf is its 
ct-component (a = x, y, z), and (ij) runs over all nearest- 
neighbor pairs. The exchange interaction Jy takes on 
either +J or — J with the same probability of 1/2. We 
assume that the anisotropy comes from pseudo-dipolar 
couplings and impose the restriction Df? — Djf = D^ a . 

We further assume that D°f are uniform random val- 
ues between —D and D. We note that the role of the 
anisotropy is to break the rotational symmetry of the 
model and to stabilize the SG order. 

We performed a MC simulation of the two-replica 
systems pfi {Si} and {Si} using an exchange MC 
algorithmcJ. We calculated the order-parameter prob- 
ability distribution Pl(q) of 



P L (q) = [(S(q-Q))} 



(2) 



where (• • •) and [• • •] mean the thermal average and the 
bond distribution average, respectively. Here Q is the 
spin overlap defined by 



a, (3 



(3) 



with q a/B = ± J2Zi S "$i- Usin S P l(q), w e obtained the 



SG susceptibility xsg and the Binder parameter g(L,T) 
which are defined by 



XSG = 3N[(q 2 }}, 
^,T)-i(ll-9g), 



(4) 
(5) 



where [(q n )] = J q n PL(q)dq. To examine whether the SG 
order occurs or not, we also calculated quantities A(L, T) 
and G(L,T) that measure the order-parameter fluctua- 
tions (OPF): 



[(a 2 ) 2 ! - [(o 2 )l 2 
A(£ ' T)= [tf : 

K<z 2 } 2 ] K? 2 )] 2 



G(L,T) 



[<<Z 4 >] - 1(<1 2 W 



(6) 
(7) 



Each of these quantities will exhibits-ajCjossing behavior 
at Tsg when the SG transition occurstJO. We performed 
the MC simulation of the model (1) with various values 
of D = 0,0.1J,0.2J,0.5J and 1.0 J. The linear sizes of 
the lattice studied here are L = 3 ~ 9. Equilibration was 
checked by monitoring the stability of the results against 



at least two-times longer runs. The numbers of the sam- 
ples were 480 for L = 3, 288 for L = 5, 96 for L = 7, and 
48 for L = 9. 
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FIG. 1. (a) Temperature dependences of the spin-glass susceptibility \ 
sizes of the lattice. Open symbols represent those for D = and closed or 
plot of xsg for D = 0. 



In Fig. 1(a), we show results of xsg m D = and 
D = 0.5J. In both the cases, xsg f° r larger L increases 
rapidly as the temperature is decreased. The same was 
true for different values of D. The finite size scaling anal- 
ysis in each D suggested the divergence of xsg for L — > oo 
at a finite, non-zero temperature. An example of the scal- 
ing plot in the case of D — is presented in Fig. 1(b). 
Hereafter, we tentatively call this temperature the SG 
transition temperature and denote Tsg- Note that Tsg 
increases with increasing D. These results are compat- 
ible with a common belief that the SG order occurs in 
d = 4 even for D ^O^and that the SG order is stabilized 
by the anisotropyB'El 

Now we show g(L, T) in different -D's in Figs. 2(a) 
- 2(c), When D = 0, g(L,T)'s for L = 3 and 5 come 
close to each other near Tsg- But they do not cross, be- 
cause the increment of g(L, T) for L = 5 is suppressed 
below Tsg- This suppression is not released for larger 
L. Therefore we believe that, in contrary to the previ- 
ous suggestiontj, g(L,T)'s for large L do not jCjoss but 
converge on some non-zero value g at T = Tsg-E3- When 
D =/= 0, the suppression is enhanced more. In particu- 
lar, for large D, g(L, T) exhibits a dip near Tsg which 
deepens as L is increased. This result is also incompat- 
ible with our naive expectation that, as D is increased, 
g(L,Tys would tend to cross with some positive g, be- 
cause the anisotropy will stabilize the SG order. 
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FIG. 2. Temperature dependences of the Binder parameter g(L,T) of 
nitude of the anisotropy D and for different sizes of the lattice; (a) D = 
and D = 1.0 J (closed symbols). Arrows indicate the SG transition tem] 
SG susceptibility. 
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We next examined G(L, T) which would exhibit a 
crossing behavior at Tsg epep-,when g{L,T) did not ex- 
hibit the crossing behaviorE2E3. Here, in Figs. 3(a) and 
3(b), we show G(L,T)'s for different L in the cases of 
D = and D = 0.5J, respectively. When D = (and 
also D < 0.2 J), G(T,T)'s for large L(> 5) seem to come 
together near Tsg- This property becomes more promi- 
nent in large D, and the data for L = 3 joins. A sim- 
ilar crossing has also been found in the other quantity 
A(L,T). 



On the other hand, in the case of small D, no dis- 
tinct dip but a bending of g(L,T) was seen around Tsg- 
We may also explain this behavior based on a plausible 
assumption that 5(00, Tg G ) takes a negative value or a 
small positive value due to the occurrence of the one- 
step RSB-like state. When D — 0, the rotational sym- 
metry recovers and g(oo,T) might exhibit some different 
property. We think, however, that g(oo, T^ G ) also takes 



a small positive value (or a negative value), because our 
data of Figs. 2(a) - 2(c) imply the presence of no gap 
between the cases of D = and D 7^ 0. Of course, we 
could not rule out the possibility that there exists some 
threshold D*, including D* = 0, below which the nature 
of 3(00, T) changes qualitatively. 

In summary, we reexamined the spin-glass phase tran- 
sition of the ±J Heisenberg model in four dimensions 
(d — 4) and gave a confirmation that the SG transition 
really occurs even when the anisotropy is absent D = 0. 
FIG. 3. Temperature dependences of the order-parameter fluctH#S&*¥^MX^ 
for different sizes of the lattice and for different magnitude of th^ar^$flto>rfit^ T (fl')t : }$ ff«Qiflqd r r(liOiT& ^r&5\/inACTOM^iwfe^ifo; 
the SG transition temperature Tsg estimated from the scaling plo^ r gf^h^S(^ (S^sqejiiqbijll^f. f rorn the crossing of G(L TVs 

(and also A(T,T)'s), as well as from the divergence of 



EPS File FIG3N.eps.epsf not 
found 



We have calculated four quantities xsg, 9(L,T), 
G(L,T) and A(L,T) to examine the SG phase transi- 
tion. All the quantities except for g(L, T) suggested the 
occurrence of the SG order below Tsg- Considering the 
facts that G(L,T) and A(L,T) can be used to deter- 
mine the value of Tsg of such SG models in which Tsg 
is hardly, jietermined by the usual crossing behavior of 
g(L, T)E3't3, and that the T-dependence of xsg clearly 
suggests the divergence of the spin correlation length at 
Tsg, it is natural to conclude that the SG order occurs 
below Tsg- Therefore, the absence of the crossing of 
g(L,T)'s for finite L would say nothing about the SG 
phase transition in this model. 

What occurs in g(L, T) at Tsg for L — > 00? The occur- 
rence of a dip of g(L, T) near Tsg for large D (D = 0.5 J 
and 1.0 J) is suggestive. The dip deepens as L is increased 
implying a negative divergence of g(oo, TJq). It is known 
that, in the infinite-range p-state Potts glass, the value of 
3(00, Tp G ) changes continuously from 1 for p — 2 (Ising 
model) to — 00 for p — . 4. due to the occurrence of the 
one- step RSB for p > 2BB. Therefore we suggest that 
3(00, T) takes —00 or some very large negative value at 
Tg G due to the occurrence of a one-step RSB-like state. 
In fact, as shown in Fig. 4, the order parameter distribu- 
tion function Pl(q) for D = 1.0 J exhibits a double peak 
at low temperatures^! 



the SG susceptibility xsg- This fact was quite interest- 
ing, because it threw a doubt on the common belief that 
the most reliable evidence of the occurrence of the SG 
order is the crossing of g(T,T)'s. We speculated that the 
low temperature phase of the vector SG model is char- 
acterized by one-step RSB-like state and the crossing of 
g(L, T)'s for finite L is absent. 
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The crossing behavior of g(L,T) at T = Tsg comes from 
the combination of (A) the scaling hypothesis and (B) the 
assuption g(oo,T) — for T > Tsg and g(pa,T) = 1 for 
T < Tsg- Therefore when the SG phase transition occurs, 
even if the assumption (B) is not satisfied, g(L, T) will con- 
verge on some non-zero value at T = Tsg- We think a large 
finite size effect masks this property. 

' Note that a double peak was found in the chirality orde* 
parameter distribution (q) of the Heisenbergu and XYu 
SG models with D = in d — 3 in which the one-step RSB 
was suggested to occur. 
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